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8.1. Contexto 

Hemos mencionado que el problema fundamental del álgebra lineal consiste en resolver un sistema de 
ecuaciones lineales Ax = b y hemos introducido varios conceptos sofisticados para manejar la solución y el 
análisis. El concepto de espacio generado resuelve el problema de la consistencia: El sistema es consistente si 
y sólo si b G C( A) (el espacio generado por las columnas de A). Otro concepto importante es el concepto 
de dependencia lineal y resuelve el problema de la unicidad: Si el sistema A x = b es consistente, entonces 
Ax = b tiene infinitas soluciones si y sólo si las columnas de A forman un conjunto linealmente dependiente. 
Los conceptos de dimensión y de base de un espacio lineal permiten tener una idea qué tan grande es el conjunto 
solución y cómo generarlo: Si el sistema A x = b es consistente, entonces la fórmula de todas las soluciones es 
y = y D + ns(A) donde y G es una solución particular y ns(A) es el espacio nulo de A, es decir, el conjunto de 
todos los vectores que satisfacen Ax = 0. De hecho, si SS = {xi, . . . ,x r } es una base para ns(A), entonces 
la solución general será 

r 

y = y o + p • 

i=l 

y esta será la solución general en su forma más reducida posible. El caso consistente ha dado origen a toda la 
teoría que hemos visto hasta ahora. Pasemos ahora a la inconsistencia. De alguna forma debemos cambiar la 
pregunta por que el tema se agotó: Si es Ax = b inconsistente, entonces no hay solución. Una forma adecuada 
de reformular la pregunta es: No habiendo solución para Ax = b, ¿qué es lo más cerca posible que podemos 


estar de una solución ? Esto involucra dos conceptos hermanados: distancia y perpendicularidad. Y ellos tienen 
como origen el concepto de producto interno o producto punto que es el tema de esta lectura. 

8.2. Introducción 

En esta lectura veremos cómo definir un producto punto o producto interno en espacios de matrices (pues 
tenga presente que nos interesa el paso de Ax = b a AX = B también en el caso inconsistente) . Veremos 
que aunque en un principio la definición del producto punto apartir del uso de la traza sea extraña , esta 
definición es de lo más cómoda porque hace coincidir el producto punto entre matrices con el producto punto 
tradicional entre vectores cuando las matrices son vectorizadas. La estructura clásica de la construcción del 
Algebra Lineal referente al producto interno es: definir un producto interno, probar ciertas propiedades básicas, 
estas propiedades incluyen en concepto de ortogonalidad, y entonces mostrar que se puede definir una norma a 
partir del producto interno; para llegar la definir una distancia con la norma se requiere probar la desigualdad 
de Cauchy-Schwarz. Cuando se dice que se tiene una distancia en el espacio lineal se tiene el concepto de 
cercanía y el de error. Las propiedades importantes que deben ser deducidas para que la distancia definida 
sea cabalmente una distancia son de identidad (que si la distancia entre dos objetos es cero, entonces los 
objetos son idénticos), de simetría (que la distancia medida de un objeto a otro es independiente de cual 
de los dos sea el punto de referencia), la desigualdad del triángulo (la distancia entre dos puntos no excede 
la suma de las distancias de uno de ellos a otro intermedio mas la distancia de ese intermedio al segundo 
punto; y habrá igualdad cuando el punto intermedio esté en el segmento que une ambos puntos). El proceso 
de Gram-Schmidt será la clave para encontrar una base para C(A) desde donde será fácil encontrar el punto 
más próximo de C( A) a b y allí a la mejor solución posible a Ax = b. 

8.3. Un producto punto mediante la traza 

El producto punto entre matrices se define apartir del producto, la transpuesta y la traza entre y de 
matrices. En apariencia complejo, veremos su conveniencia y facilidad posteriormente. 

Definición 

Sean A y B dos matrices rn x n, el producto interno o producto punto de A y B se define como: 


A • B = tr (A' • B) 


(1) 


Ejemplifiquemos el cálculo de productos punto. 

Ejemplo 1 

Determine A • B si 



' 3 

-1 ' 


' 1 

-1 ' 

A = 

2 

0 

, B = 

0 

1 


5 

-2 


4 

-3 


Solución 


A' B = 




' 1 

-1 ' 




3 2 

5 


= 

23 

-16 ' 

0 

1 

-1 0 

-2 



-9 

7 




4 

-3 





Por tanto: 


A • B = tr 


23 -16 
-9 7 


= 23 + 7 = 30 


Toca el turno al alumno para ejercitar el concepto del producto punto entre matrices de la misma dimensión. 

Ejercicio 1 


2 



Para las matrices siguientes calcule A • B. 


1. A = 


2 -1 
3 1 


y b 


1 o 

-1 o 


2. A = 


2-13 
1 O 1 


yB 


i o -i 
O 2 1 



' 2 

1 

i-H 

1 


i 

i 

O 

3. A = 

3 

1 

y b = 

-i 

0 


1 

O 

1 _ 


i 

LO 

1 _ 


8.4. Propiedades del producto punto 

El camino tradicional para ver que un producto es efectivamente candidato a producto interno con escalares 
reales es probar que es: conmutativo, no negativo, y bilineal: 


Teorema 8.1 

Sean A, B y C matrices mxny sea k un escalar. Entonces el producto punto cumple las siguientes 
propiedades: 

1. Simetría: A • B = B • A 

2. No negatividad: A • A > 0, habiendo igualdad ssi A = 0 

3. Proporcionalidad: ( kA ) • B = k (A • B) 

4. Distributividad: (A + B)»C = A*C + B*C 

5. 0* A = 0 


Demostración 


1. Directo de la definición: 

A • B = tr(A'B) 


= tr((A'B/) 

= tr(B'A) 


= B • A 

2. De la sección anterior: 

A. A = tr(A'A) 


= £í=i E j= i % > o 

Además, A • A = 0 ssi 

m n 




i = 1 i 1 


ssi a,ij = 0 para todo i = l...myj = l...nssiA = OH 
No sea impaciente; Ud. participará en el festín de mini-demostraciones: 

Ejercicio 2 

Demuestre el punto 3 del teorema anterior: 

Sean A y B matrices mxny sea k un escalar. Pruebe que: 

(kA) • B = k (A • B) 


Sugerencia 

Utilice la definición de •, la propiedad de la transpuesta, y la propiedad 1 del Lema 5.1. 
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Ejercicio 3 


Demuestre el punto 4 del teorema anterior: 

Sean A, B y C matrices m x n. Pruebe que: 

(A + B).C = A*C + B*C 


Sugerencia 

Utilice la definición de •, la propiedad de la transpuesta, y la propiedad 2 del Lema 5.1. 


Ejercicio 4 

Demuestre el punto 5 del teorema anterior: 

Sea A una matriz m x n. Pruebe que: 

0* A = 0 


8.5. La norma de una matriz 

Habiendo probado las propiedades mencionadas del producto interno nuevo, el siguiente paso es definir la 
norma de una matriz. En el espacio de matrices hay un buen número de normas igualmente útiles, la norma 
que se define aquí es la norma-2 de una matriz. 

Definición 

La norma de la matriz A m x n se define como: 


A|| 2 = || A|| = Va • A 


(2) 


Observe que por el punto 2 del teorema anterior, A • A > 0; así que la raíz cuadrada entregará un real mayor 
o igual a cero. Hagamos un ejemplo del cálculo de la norma de una matriz. 

Ejemplo 2 

Determine la norma de la matriz: 


Solución Como 


A'A 


3 

-1 



Así 


A • A = tr 


38 -13 
-13 5 


= 43 


y por tanto, ||A|| = \/43 ■ 

Es su turno de practicar el cálculo de la norma de una matriz. 

Ejercicio 5 


4 



Para las matrices siguientes calcule su norma. 


1. A = 


2 -1 
3 1 


2. A = 


2-13 
1 0 1 


3. A = 


2 -1 
3 1 

0 1 


Un hecho importante sobre la norma de una matriz es: 


1/2 


| A|| = tr (A' A) = [ 

*=1 3 = 1 


Lo cual hace coincidir la definición tradicional de norma de un vector con la norma de 
matriz es utilizada para formar un solo vector formado por los renglones de ella. Algo 
observar es que el producto punto entre matrices usando la traza equivale 
tradicional aplicado a los vectores obtenidos de vectorizar las matrices. 
Ejemplo 3 

Verifiquemos la afirmación en el caso de matrices 2x2. 

Comprobación 

Digamos que 


an 

ai2 

B = 

' 611 

¿>12 

021 

022 

^21 

t>22 


una matriz cuando la 

más importante de 
al producto punto 


Así 


A'B 


Oll&ll + 021^21 



012^12 + O 22^22 


De donde 


A • B = (an&ii + 021621 ) + (012612 + 022622 ) 


8.6. Propiedades de la norma de matrices 

Habiendo definido la norma de una matriz, lo siguiente es probar algunas propiedades básicas. 


Teorema 8.2 

Sea A una matriz m x n y k un escalar: 

1. 1 1 A 1 1 > 0, además ||A|| = 0 ssi A = 0 

2. || - A || = || A || 

3. \\kA\\ = \k\ || A || 


Demostración 

1. De lo visto 


A 1 1 = V A» A 


\ 1/2 

m n \ 

¿¿M 2 >0 

i= 1 .7 = 1 / 


5 


Y también de lo mismo se deduce que ||A|| = 0 ssi A = 0. 

3. Por las propiedades conocidas: 

||fcA|| = \¡ (kA) • (kA) = Vk 2 A • A = \k\ V A • A = \k\ ||A 

El segundo punto es el más sencillo y le queda hacerlo como ejercicio ■ 

Ejercicio 6 

Demuestre el punto 2. del teorema anterior: 

Sea A una matriz m x n. Pruebe que: 


- Alt = II Al 


Sugerencia 

Use directamente la definición de la magnitud de una matriz. 

8.7. Desigualdad de Schwarz 

En la línea tradicional, el siguiente paso para ver que nuestro producto punto llega a definir una distancia 
a través de la definición de la norma, es probar que nuestro producto punto y la norma cumplen una relación 
de desigualdad llamada la desigualdad de Cauchy- Schwarz. 

Teorema 8.3 

Sean A y B matrices m x n: 

( 3 ) 

Además, la igualdad se tiene si y sólo si una matriz es un múltiplo de la otra. 

Demostración 

En esta demostración es importante que vigile el rol de los factores en los productos: en los productos • los 
factores en ambos lados deben ser matrices, mientras en los productos • el factor a la izquierda es escalar y el 
factor a la derecha debe ser una matriz. Para cualesquiera escalares x y y: 

0 < (x ■ A — y ■ B) • (x ■ A — y ■ B) 

= x 2 || A|| 2 — 2xy (A • B) + ?/ 2 ||B|| 2 

Suponiendo B / 0, tomamos x = ||B|| y y = (A • B)/||B|| en la fórmula anterior y simplificando muchos 
términos obtenemos: 

0 < 1 1 B 1 1 2 || A|| 2 - (A • B) 2 

De donde se obtiene: 

(A • B) 2 < ||B|| 2 ||A|| 2 

y tomando raíz cuadrada se obtiene la conclusión. Recuerde que un real k, Vk 2 = |fe| y no que Vk 2 = k pues 
para los valores de k negativos el signo se pierde al elevar al cuadrado. Además, si A = k ■ B se tiene que 
|| A || = \k\ ||B || y así 

|A • B| = |jfe| (B • B) = \k\ || B || 2 = 1 1 A 1 1 || B || 

probado que la desigualdad se cumple con igualdad en este caso. Por otro lado, si|A»B| = ||A||||B|| definimos 

C = ||B|| 2 -A-(A*B)-B 


A • B| < ||A|| • ||B 


6 



se puede confirmar que desarrollando el producto término a término: C • C = 0, implicando que C = 0, de 
donde 1 1 B 1 1 2 • A = (A • B) • B si suponemos B / 0 deducimos que 


(A.B) 

1 1 B 1 1 2 


B 


lo que dice que A es un múltiplo de B. Si B = 0, claramente B = 0 • A y la desigualdad se cumple con 
igualdad ■ 

Ejercicio 7 


Sean A y B dos matrices m x n. Demuestre que 


A + Bll < IIAII + II B 


y que existe igualdad sólo cuando una matriz es un múltiplo escalar de la otra. 

Sugerencia 

Utilice 

|| A + B|| 2 = (A + B) • (A + B) 

Desarrolle tal producto punto. Después utilice la desigualdad de Schwarz. Finalmente, tome raíz 
cuadrada. Para la igualdad se requiere igualdad en la desigualdad de Schwarz. 


Lema 8.4 

Sean X y Y dos matrices m x n. Entonces X y Y son ortogonales, es decir, cumplen X • Y = 0 si 
y sólo si 

||X± Y|| 2 = ||X|| 2 + || Y|| 2 


Demostración 

Tenemos en general que 

1 1 X ± Y 1 1 2 = (X ± Y) • (X ± Y) 

= X.XÍ2X.Y+Y.Y 
= ||X|| 2 + ||Y|| 2 ±2X*Y 

Por tanto, si X • Y = 0 entonces de la fórmula anterior se deduce la igualdad que queremos probar. Por otro 
lado, si la fórmula se cumple nuestra fórmula indica que X • Y = 0. Es decir X y Y son perpendiculares. ■ 


8.8. La distancia entre dos matrices 


Habiendo probado la desigualdad de Cauchy-Schwarz se puede construir una función distancia entre el 
espacio de matrices: 

Definición 

Sean A y B dos matrices m x tí la distancia entre matrices A y B se define como: 


6 (A, B) = ||A-B 


Realicemos un ejemplo simple para ilustrar este concepto. 
Determine la distancia entre las matrices: 


' 3 

-1 ' 


' 1 

-1 ' 

2 

0 

B = 

0 

1 

5 

-2 


1 

-2 


(4) 
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Solución Si 


Así 


C'C 


c = A-B = 


2 0 
2 -1 
4 0 


2 2 
0 -1 


4 

0 


2 0 
2 -1 
4 0 


24 -2 
-2 1 


Í(A,B) = ||C 


VcTc 


^tr(C'C) = \/24+T = 5 0 


Ejercicio 8 

Determine la distancia entre las matrices: 



1 

co 

1 

I— 1 

1 


l 

l ' 

1. A = 

2 

1 

II 

ffl 

-l 

l 


_ 5 

1 

CN 


l 

to 

l 


2. A = 


3 -1 
5 2 


yB 


-i i 
1 2 


Ejercicio 9 

Sean A y B dos matrices de iguales dimensiones. Pruebe que: 

Í(A,B) = Í(B, A) 


Sugerencia 

Use directamente la definición de distancia entre matrices. Utilice también la propiedad conmutativa 
del producto punto. 

Ejercicio 10 

Sean A y B dos matrices de iguales dimensiones. Pruebe que: 

Ó(A, B) > 0 

Además, pruebe que hay igualdad si y sólo si A = B. 

Sugerencia 

Utilice la definición de distancia entre matrices y la propiedad 1 del teorema 6.2 sobre la norma 
para matrices. 


Ejercicio 11 

Sean A, B y C dos matrices de iguales dimensiones. Pruebe que: 

¿(A, B) < 5(A, C) + <5(C, B) 

Vea que existe igualdad si sólo si existe un escalar t (0 < t < 1) tal que C = t A + (1 — t) B. Esta 
desigualdad se llama Desigualdad del triángulo. 

Sugerencia 

Utilice como válido el resultado del ejercicio 7 con A como A — C y con B como C — B. 
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Ejercicio 12 

Sean A, B y C dos matrices de iguales dimensiones. Pruebe que: 

Ó(A, B) = Ó(A + C, B + C) 

Esta propiedad indica que la distancia entre matrices se preserva ante traslaciones. 

Sugerencia 

Directamente de la definición de distancia. 


8.9. Ángulo entre matrices 
Definición 

Sean A y B dos matrices m x n, no nulas, el ángulo entre matrices A y B, 9, se define como 


Aquí 9 debe cumplir: 0 < 6 < n. 


eos (9) 


A • B 

Al II B 


Ejemplo 4 

Determine el ángulo en radianes entre las matrices 


Solución Como 



A • B = —2, || A 


a/43, ||B 


V22 


Entonces 

_2 

cos(0) = , . « —0.065, por tanto 9 ss 1.63 radianes<> 

^ 43^22 


(5) 


Ejercicio 13 

Determine el ángulo en radianes entre las matrices 



Ejercicio 14 

Determine el valor de x para que el ángulo en radianes entre las matrices 


A = 


3 

5 


-1 

-2 


B = 


-4 -1 

1 x 


sea de 0.1. 

Sugerencia 

Entróle sin miedo: A • B es una expresión lineal en x. |B¡ es una raíz cuadrada de un polinomio 
cuadrático en x. Cuando se iguala: 

. . A • B 

cos(O.lrad) = 

A || B || 

después de elevar al cuadrado y multiplicar por denominadores queda una ecuación cuadrática 
en x. 
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8.10. Ortogonalidad y espacios generados 

Sean A y B dos matrices m x n, se dice que son matrices ortogonales si se cumple 

A • B = 0 


Un conjunto {Ai, A 2 , . . . , A*,} de matrices m x n se dice conjunto ortogonal si A ¿ • Aj = 0 para i 7 ^ j El 
conjunto anterior se dice ortonormal si es ortogonal y además cada matriz tiene norma 1. Una vez que se 
tiene un conjunto ortogonal de elementos no cero 


es posible convertirlo a uno ortonormal 


definiendo: 


{Ai, A 2 , . . . , A*,} 
{Bi, B 2 , . . . , Bfc} 

B i = -¡j-T— ñ A* 

II -^-2 || 


Ejercicio 15 

Determine el valor de x para que las matrices sean ortogonales. 


A = 


3 

5 


-1 

-2 


B = 


-4 -1 

1 x 


Sugerencia 

Al hacer A • B = 0 queda una ecuación lineal en x. 


Ejercicio 16 

Determine qué condición deben cumplir las variables para que A A B. 


A = 


3 

5 


-1 

-2 


B = 


x y 
z w 


Sugerencia 

Al hacer A • B = 0 queda una ecuación lineal. 


Ejercicio 17 

Sea A una matriz m X n muestre que el conjunto 

A^ = {X matriz m x n|A A X} 


es un espacio lineal. 

Sugerencia 

Al hacer A • X = 0 queda una ecuación homogénea. Las soluciones a los sistemas homogéneos son 
espacios lineales. 


Ejercicio 18 
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Determine el espacio perpendicular a la matriz: 


1 1 

-1 1 


Escíbalo como un espacio generado. 

Sugerencia 

Tome una matriz X 2 x 2 con incógnitas. Al hacer A • X = 0 queda una ecuación homogénea con 
4 incógnitas; una de ellas es fija y las otras son libres. Despeje la fija y póngala en función de las 
libres. En lugar de la notación de vector en la solución general, prefiera la notación de matriz. 


Ejercicio 19 

Sea V un subespacio lineal de matrices m x n muestre que el conjunto 

V ± = {X matriz m x n|Y _L X para toda Y E V} 


es un espacio lineal. 

Sugerencia 
Muestre que 

■ no es vacío: la matriz de ceros está (compruébelo), 

■ es cerrado bajo la suma (compruébelo), y 

■ es cerrado bajo el producto por escalares (compruébelo). 


Ejercicio 20 


Determine el espacio perpendicular al espacio generado por las matrices matrices: 


A = 


1 1 

-1 1 


y b 


i -i 
i i 


Sugerencia 

Tome una matriz X 2 X 2 de incógnitas: 

■ A • X = 0 da una ecuación en las incógnitas 

■ B • X = 0 da otra ecuación en las incógnitas 

Resuelva el sistema homogéneo, encuentre la solución general y déle la notación matricial. 
Veamos ahora un par de resultados teóricos: 

Lema 8.5 

Un conjunto ortogonal de matrices no nulas es linealmente independiente 


Demostración 

Si {Ai A 2 , . . . , Ai,} es un conjunto ortogonal de matrices no nulas y si se cumple: 


c\ Ai + C2A2 + • • • + CfcAfc — 0 
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veamos que los coeficientes c¿ son cero. Haciendo el producto punto con A¿: 

A¿ • (ci Ai + c 2 A 2 H b c k A fc ) = A¿ *0 = 0 

Desarrollando el primer miembro: 

Aj • (ciAi) + A¿ • (c 2 A 2 ) H f A ¿ • (c fe A fc ) = 0 

Como el conjunto es ortogonal entonces A¿ • A j = 0 para i ^ j, la expresión anterior queda: 

c¿A,¿ • A¿ = 0 

Puesto que ninguna matriz A ¿ es nula, A ¿ • A¿ / o ,1a expresión anterior queda: 

c¿ = 0 

Así, cada coeficiente en la combinación lineal es cero. Por tanto, el conjunto de matrices es linealmente 
independiente. ■ 

Cuando se quiere verificar que una matriz es ortogonal al generado por un conjunto de matrices enfrentámos 
un dilema singular, cómo hacer pues en el generado resulta en conjuntos infinitos? El siguiente resultado es 
importante porque reduce una verificación infinita a una verificación sobre el conjunto generador: 

Lema 8.6 

Sea A una matriz y V = Gen{Bi, .... B^} en M mxn . Entonces: A es ortogonal a todo V si y sólo 
si A es ortogonal a cada B¿. 


Demostración 

Suficiencia : Si suponemos que A es ortogonal a todo V y sabiendo que cada B¿ está en V, entonces es inmediato 
que A es ortogonal a B¿. 

Necesidad : Supongamos ahora que A es ortogonal a todo B¿, es decir que A • B¿ = 0 para todo i = 1, . . . , k. 
Veamos que entonces es ortogonal a cualquier X de V. Tomemos un X cualquiera de V. Como V está generado 
por los B¿, entonces deben existir escalares c¿ tales que 


Por tanto, 


Por tanto, A es ortogonal a X ■ 


X = ^ c fc • B¿ 


A • X = A.(Ef=iC fe -B ¿ ) 
= E¿= i A • (cfc • B¿) 
= E¿=i C ¿(A • B¿) 

= £¡U«(0) 

= o 


8.11. Bases Ortonormales 

El siguiente resultado es importante porque permite cambiar la base de un espacio generado por otra que 
es ortogonal. La misma demostración del resultado da el algoritmo de conversión. 


Teorema 8.7 
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Sea un conjunto {Ai, A 2 , . . . , A*,} de matrices m x n linealmente independiente. Entonces existen 
escalares únicos < j) tales que: 


Bi 

= A, 




b 2 

= a 2 

~ Zl,2Bi 



B i 

= A-j 

- xijBi 

- x 2 ,jB 2 


B k 

= A k 

- íCgfcBi 

- X 2 ^2 ~ ■ 

x k— k — 1 


forman un conjunto ortogonal. 


Demostración 

Definimos los valores escalares xjj como precisamente se requieren: 


B; 


x i,j — 


B, • B, 


(6) 


Ud. puede recordarlos si ubica que el primer subíndice de x t j se relaciona a las matrices por construir B¿ y que 
ésta es la que aparece en el denominador con el producto punto con si misma. Un punto importante a observar 
es que los {Ai, A 2 , . . . , A j } son combinaciones lineales de los {Bi, B 2 , . . . , Bj} y por consiguiente el conjunto 
de los B.s es linealmente independiente. Para verificar que el conjunto es ahora ortogonal hay que realizar 
los productos punto. Conviene hacer una demostración inductiva de la afirmación. La afirmación es que el 
conjunto {Bi, . . . , B&} es ortogonal. Esta afirmación sería la afirmación última de la cadena de afirmaciones 
{Bi, ■ • • , B.¿} para i = 2, 3, . . . , k. La base inductiva para i = 2 se comprueba revisando que Bi y Bo son 
ortogonales: 

A 2 • Bi 


B2 • Bi = A2 • Bi — xioBi • Bi = A2 • Bi — 


-B i • B i =0 


Bi • Bi 

Supongamos que {Bi,...,B t } es un conjunto ortogonal para 1 < i < k — 1, veamos {Bi, . . . , B¿, B¿+i} 
también lo será. Para ello basta probar que B¿ + i es ortogonal a cada Bj para 1 < j < i. Como 


B¿+i — A¿_|_i ^ ^ x^j-i-iBj 


1=1 


Así 


B¿+i • B j = A¿+i - J2i = 1 x i,i+iRi 


1 B, 


= A¿ + i • Bj - ^]=i x i,i+i&i • Bj 


= A 

= 0 


í ■ 1 • B j Xj.i ■ 1 1 {y • Bj 


Note que en el penúltimo paso fue requerida la base inductiva: que Bj era ortogonal a sus compañeros en 

{B 1; . . . , B.¿}.« 

Ejemplo 5 

Aplique el algoritmo anterior al conjunto formado por las matrices: 


Ai = 


1 0 
1 0 


• A 2 = 


1 -1 

1 o 


, A;¡ = 


1 1 
0 2 


Solución 


B, = Ai 
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' f1 7*jT 






■ El □ 1 O] + al 

■El - 1 1 0] ■» a2 

■ El 1 B 2] -> a5 


[10 10 ] 
[1 “I 1 O] 


LINEAL 


DEG EXflCT FLNC 1/ í 


Figura 1 : Captura de la matrices vectorizadas del ejemplo 7. 


Hjmjr í FSt T FHt y FE Y FSt V 

■|Hlgebra|Calo|0ther|Pr9nl0|Clean Up| 


- 1 1 0] ■» a2 

■ El 1 0 2] a3 

■ al + bl 

. dotP(a2 , bl) 


■ a3 - 


dotPíbl , bl) 
dotP(a3, bl) 
dotPCbl , bl) 


[1 ‘I 1 O] 


[1 

El 


bl+b2 [0 -1 

dotP(a3, b2) 


bl 


dotP(b2 , b2) 
[ 1/2 0 - 


0 0 ] 

■ b2 * 
1/2 2 ] 


LINEAL 


DES EmIÍCT FLNC B. í'> 


Figura 2: Obtención del conjunto ortogonal del ejemplo 7. 


Realizando los cálculos: 


Bi • Bi = 2 y A 2 • Bi = 2; X 12 = 1 
B 2 = A 2 — xi 2 Bi = 


B 2 • B 2 = 1; A 3 • Bi = 1; A 3 • B 2 = — 1; 


Obtenemos: 


£13 = 1 / 2 , x 23 = - 1 ; 


Y así: 


B 3 = A 3 — X 13 B 1 — x 23 B 2 


1/2 0 

- 1/2 2 


Estos cálculos pueden hacerse en la calculadora TI, para ello habrá que vectorizar las matrices y seguir las 
fórmulas: en la figura [T] aparecen capturadas las matrices como vectores y en la figura [2] aparecen los cálculos 
del algoritmo de ortogonalización. 

Ejemplo 6 

Aplique el algoritmo anterior al conjunto formado por las matrices: 



' 1 

0 ' 


' 2 

0 ' 


■ 1 

1 ' 

Ai = 

, A 2 — 


, A 3 = 

1 

0 


2 

0 


0 

2 


Solución 


B 1 = A 1 


Bi • Bi = 2 y A 2 • Bi = 4; x 12 = 2 
B 2 = A 2 -xi 2 Bi= ° ¡j 

este cálculo indica que vector A 2 es combinación lineal de los vectores anteriores a él y por tanto el conjunto 
original es linealmente dependiente. Debemos descartar al vector A 2 . Realizando solo el cálculo: 


A3 • Bi = 1; 
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Obtenemos: 


*13 = 1/2 


Y así: 


B3 = A3 — .T13B1 


1/2 1 
- 1/2 2 


Ejercicio 21 


Aplique el algoritmo anterior al conjunto formado por las matrices: 


Ai = 


1 0 
1 1 


5 A 2 


1 -1 

1 1 


■ A 3 


1 o 
0 2 


Sugerencia 

Siga el algoritmo; no hay de otra. 

8.12. Existencia de Bases Ortonormales 

EL resultado anterior se conoce como el procedimiento de Gram-Schmidt, y permite cambiar una base por 
una base ortogonal, de donde es fácil obtener una base ortonormal: 

Teorema 8.8 

Todo espacio lineal de matrices m x n posee una base ortonormal 


Demostración 

Sabemos que todo espacio lineal posee una base, digamos {Ai, . . . , A¿.}. Sean lo Bj los las matrices obtenidas 
de aplicar el proceso de Gram-Schmidt. Los espacios generado por los B¿ y los A ¿ son iguales pues las matrices 
B, son combinaciones lineales de los vectores A ¿ y viceversa. Así, los conjuntos generados son los mismos. Si 
las matrices A¿ forman un conjunto linealmente independiente, también lo debe ser el conjunto formado por 
los Bj. Por tanto, el conjunto formado por los B¿ es también una base para el espacio generado por los A ¿. 0 


Una ventaja inmediata de una base ortonormal {Ai, A 2 , ... . , A*.} es que si A es una combinación lineal 
de las matrices anteriores entonces: 

A = (A • Ai) Ai + (A • A 2 ) A 2 H h (A • A fc ) A k 

es decir, que no hace falta utilizar el proceso de eliminación gaussiana para determinar los coeficientes, si no 
simplemente hacer un producto punto con el vector correspondiente. 

Una consecuencia también inmediata es que se puede completar un conjunto linealmente independiente y orto- 
gonal a una base ortogonal en un espacio lineal de matrices mxn. Una forma conviente es aplicar el proceso de 
Gram-Schmidt al conjunto li y ortogonal aumentado con una base cualquiera del espacio y detener el proceso 
hasta obtener el primer vector cero. El resultado esperado debe ser el conjunto inicial adicionado con algunos 
vectores. Este conjunto debe ser una base para el espacio completo. 
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8.13. Proyección ortogonal 


Teorema 8.9 

Sea Y una matriz m x n y un espacio lineal Y de dimensión r, ambos dentro de un espacio lineal 
U. Entonces, existe una única matriz Z en Y tal que (Y — Z) _L Y. Si r = 0 entonces Z = 0, y 
si r > 0 entonces Z se puede expresar como 

Z = ciXi T • ■ • T c r X r , 

donde {Xi, . . . , X r } forman una base ortonormal de Y y c¿ = Y • X,¿ para i = 1, . . . ,r. Además, 
Z = Y si y sólo si Y G Y. La matriz Z se llamará la proyección ortogonal de Y sobre Y y cumple 
que d(Z, Y) < d(X, Y) para toda X en Y y hay igualdad si y sólo si Z = X. 


Demostración 

Si r = 0 entonces dim(Y) = 0, y por tanto V = {0}. Para Z = 0 se tiene (Y — Z) _L V. Y es claramente la 
única matriz en V que cumple esto. 

Si r > 0 sea {Xi, . . . , X r } una base ortonormal de Y y definamos c¿ = Y »X,¿ para i = 1, . . . , r y Z = Ya=i 
Claramente, Z G Y y 

(y-¿ c ,x,). x í = y.x í -c ¡ = Cj - Cí =o 

para cada j = 1, . . . , r. Y por tanto, (Y — Z) _L Y. Si ahora X G Y y cumple (Y - X) 1 Y entonces: 

(X-Z).(X-Z) = (X — Y + Y — Z) • (X — Z) 

= — (Y — X) • (X — Z) + (Y — Z) • (X — Z) 

= -0 + 0 = 0 

Esto debido a que X — Z G Y y a que (Y — Z) • (X — Z) = 0. Por tanto, X — Z = 0 y de allí que X = Z, 
haciendo que Z sea el único vector en Y que cumple (Y — Z) _L Y. Si X es una matriz en Y entonces X — Z 
también está en Y y por tanto X — Z será ortogonal a Y — Z. De donde deducimos que 

||X - Y 1 1 2 = ||X - Z|| 2 + || Y - Z|| 2 

de donde 

d(X, Y) 2 = d(X, Z) 2 + d{ Y, Z) 2 

y por tanto d{ Z, Y) < d(X, Y). Y esto se cumple con igualdad si y sólo si X = Z ■ 

Ejercicio 22 

Considere el espacio lineal formado por todas las soluciones al sistema homogéneo: 

x + y + z — w = 0 

x — y — z + w = 0 


y el vector d =< 1, 3, 2, 1 >. 

■ Usando el orden primero x, luego y , luego z, y por último w , encuentre una base para el 
espacio solución. 

■ Ortogonolice la base encontrada. 
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■ Usando la base encontrada, determine la proyección ortogonal de d sobre tal espacio. 

■ Usando el orden primero y, luego w, luego y, y por último z, encuentre una base para el 
espacio solución. 

■ Ortogonolice la base nueva base. 

■ Usando la nueva base encontrada, determine la proyección ortogonal de d sobre tal espacio. 


8.14. Descomposición QR 

Sea A una matriz m x k de rango columna completo. Esto es, una matriz de rango k. Simbolicemos por 
ai, ao, . . . , afc las columnas de la matriz A. Por el resultado anterior, deben existir escalares Xi¿(i < j = 1, . . . , k) 
tales que los vectores columna bi, . . . ,b¿, definidos en forma recursiva por las igualdades: 

bi = ai 

b 2 = a 2 — aqobi 

kj = a j ~ — • • • — Xj-ijbj-i 

bfc a¿, xpfcbi ■ %k— i,fcbfc— i 

o por las igualdades: 

ai = bi 

a 2 = b 2 + aqobi 


aj = bj + xijbi + ■ ■ • + xj-ijbj-i 


— bfc + xi 5 fcbi + ■ ■ ■ Xk— i^bfc—i 

forman un conjunto ortogonal. Sea B la matriz cuya columna i es el vector b¿, y sea X la matriz k x k triangular 
superior unitaria cuyo elemento (i. j) es Xi¿ entonces: 

A = BX 


Si se define la matriz D como 

D = diag ( 1 1 b 1 1| _1 , ||b 2 || -1 , . . . , ||b fc || _1 ) , 
E = diag (||bi || , ||b 2 ||, . . . , ||b fe ||) 

y como 

Q = BD y R = EX 


entonces 

A = QR 


es una factorización de A en dos matrices: la primera formada por columnas que son ortonormales y la 
segunda triangular superior con elementos de la diagonal principal positivos. Esta factorización se conoce 
como la factorización QR de A y es muy importante en métodos numéricos. 

Ejemplo 7 

Determine la factorización QR de la matriz: 


A = 


1 1 
1 2 
1 3 
1 4 
1 5 
1 6 
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Solución : Trabajemos con las columnas de A: 


a x = (1, 1,1, 1,1,1)', a 2 = (1,2, 3, 4, 5,6)' 

Así: 

bi = ai 

Calculando: 

bi • bi = 6 y a 2 • bi = 21 

Obtenemos: 

*12 = 7/2 

Y así: 

b 2 = a 2 — si 2 bi = (-5/2, -3/2, -1/2, 1/2, 3/2, 5/2)' 

Calculando: 

b 2 • b 2 = 35/2 

Por tanto 

' 1 -5/2 
1 -3/2 

Q = [bi,b 2 ]D= ) 

1 3/2 

. 1 5/2 

’ 1/603 -1/14070 ’ 

1/603 —3/70 070 
1/603 — 1/70^/70 
1/603 1/70 070 

1/603 3/70070 

1/603 1/14070 

R _ 06 0 1 Sl2 _ 06 |03 

~ [ 0 1/2070 J [ 0 1 J “ [ 0 ^070 

Ejercicio 23 

Determine la factorización QR de la matriz: 

r 2 n 



L-i oj 

8.15. Aplicación de la factorización QR 


Teorema 8.10 

Si A es una matriz cuyas columnas son linealmente independientes y si A = QR es una factorización 
QR de A. Entonces la única solución x de Ax = b por mínimos cuadrados se expresa teóricamente 
como 

x = R _1 Q T b 


03/6 0 ' 

0 070/35 
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Ejemplo 8 

Aplique la factorización QR para resolver por mínimos cuadrados el sistema: 


" 1 

1 ' 


' 1 ' 

1 

2 


2 

1 

3 

X = 

0 

1 

4 

1 

1 

5 


1 

_ 1 

6 _ 


_ 2 _ 


Solución Por el ejemplo anterior, la factorización de la matriz 


de coeficientes es: 


A = QR = 


1/6 v/6 -l/14\/70 
l/6>/6 -3/70>/70 
l/6\/6 — l/70\/70 
1/6 v/6 l/70\/70 
l/6\/6 3/70\/70 

l/6\/6 l/14v/70 


v/6 /\/6 

0 //7Ó 


Siguiendo la estrategia descrita en el último resultado la solución por el método de mínimos cuadrados se 
calcula por la fórmula: 

x = R -1 Q T b 


Así 


q t = 


1/6V6 l/6y/6 l/6y/6 l/6y/6 1/6 v/6 l/6\/6 

-1/14V7Ó — 3/70a/7Ó -1/70 v/70 1/70^70 3/70V70 1/14V70 


R' 1 = 


l/6v/6 — l/10\/70 
0 l/35\/70 


y por tanto, 


x = R x Q T b = 


13 

15 

3 _ 
L 35 


Ejercicio 24 


Aplique la factorización QR para resolver por mínimos cuadrados el sistema: 


'21' 


' 1 ' 

-1 

1 

X = 

2 

2 

2 

0 

_ -1 

0 _ 


_ 1 _ 


8.16. Uso de la TI 

La calculadora TI 89 o TI Voyage 200 puede ser utilizada para la factorización QR. La ventaja de la rutina 
numérica programada es que da información sobre el espacio nulo de la matriz de entrada cuando las columnas 
de la matriz no son linealmente independientes. Considere por ejemplo la matriz 


Las figuras [3] y |U muestran la entrada y los cálculos en la TI. El punto a señalar es que habiendo ceros en la 
diagonal de R se indica que los vectores correspondientes en Q están en el espacio nulo de la matriz A. De 
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ú i 

■[i 

m 

■ QR a , q , r 

[f 

■ q ■ r 

Dore 

2-451 

0 

ri 2i 

.2 4 . 


LINEAL DES EKACT FUNC 5/11 


Figura 3: Descomposición QR de A en la TI. 


Ti7*U)Y fTí Y F3t y fht Y F5 Y Y 1 

■. f— |Hlgebra|Calc|other-|Pi-gri 10 Idean Up| 


1 QR a , q , r 


Done 



ri5 

2 45 ' 


Id 

□ 



2 -45 


5 

5 


2--Í5 

-45 


5 

5 . 


E 


FUNC 15/ÍO 


Figura 4: La naturaleza de los ceros en la diagonal de R. 


hecho, todos los vectores correspondientes a ceros en la diagonal de R forman la base para el espacio nulo de 
A. 

Ejemplo 9 

Determine la distancia de P(0, —2, —1) al espacio que generan los vectores: 



Con los vectores vi y V 2 formamos la matriz A a la cual le aplicamos la factorización QR: 


A = [vi V 2 ] = Q R = 


—1/7 a/14 17/1365 v/1365 
3/14 x/14 4/273 v/1365 

1/14 v/Í4 -2/105 v/1365 


v/l4 5/7 v/Í4 
0 1/7 V1365 


De donde los coeficientes de Fourier del vector que representa aP,v =< 0,-2,l > son: 


v ]q = Q v = 




y de allí que la proyección de v a C( A) es: 


161 

18 


v C (A) = Q • Mq = -# 


30 


Por tanto, la distancia de v a C( A) es 

d(v,C( A)) - 1 1 v v r(A) ;j| = 


161 


'30 


23 

390 


v/390 
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Ejemplo 10 

Determine la distancia de P(0, —2, —1, 1) al conjunto de soluciones del sistema 



/ X\ \ 

2 13 1' 

X2 

-10 2 1 

X3 

V X '4 / 


0 

0 


Este problema se resuelve como el anterior, pero debemos encontrar un conjunto generador. Para ello tenemos 
dos alternativas: o seguir el proceso descrito en lecturas anteriores o calcular lo que se conoce como el espacio 
nulo de B: Aquellos vectores x tal que Bx = 0. En general, los conjuntos generadores obtenidos por ambos 
procesos son diferentes pero ambos son base para el mismo conjunto. Si utilizamos la alternativa del espacio 
nulo tenemos: 


nullspace(B) 


Gen < 


< 

( 2 \ 
-7 

1 

i 

( 1 ^ 
-3 

0 



V 0 ) 


V i ) 



Así el problema se ha transformado en encontrar la distancia de P(= v) a C{ A) donde A es la matriz cuyas 
columnas son el generador de nullspace(B). Se deja como ejercicio comprobar que: 


\ 2 

1 1 


1/9 y/6 4/585 v/390 


-7 -3 

= QR = 

-7/18 >/6 —1/1170 \/390 


' 3 y/6 23/18 \/6 

1 

0 

1/18 y/6 -23/1170^390 


0 l/18y/m 

0 

1 


0 3/65^/390 



V C(A) 


/ 43/65 \ 
-19/65 
-47/65 
V -14/65 / 


d(v,C( A)) = l/65\/46l5 


Ejemplo 11 

Determine la matriz X que mejor resuelve el sistema AX = B para 


" 1 

2 ' 


' 1 

2 

-11 

0 

-1 

, B = 

0 

-1 

00 

1 

1 


1 

1 

to 

1 


La factorización QR de A es 


A = QR = 


1/2 y/2 1/6 x/6 

0 -1/3 y/6 

1/2 y/2 -1/6 y/6 


y/2 3/2 y/2 
0 1/2 y/6 


Por tanto, la matriz X que minimiza el error cuadrático es: 


X = r'Q'B 


1 0 10 
0 1 -19/3 


Y el error cometido con X es 25/3 ■ 
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